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Quantum-enhanced interferometry and the structure of twisted states 
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Preparation of a non-classically correlated state is the first step of any quantum-enhanced in¬ 
terferometric protocol. An efficient method is the one-axis twisting, which entangles a collection 
of initially uncorrelated particles by means of two-body interactions. Here we investigate the lim¬ 
its of the quantum improvement which can be reached with this method in realistic experimental 
conditions. We demonstrate that the usefulness of this entangling mechanism is a result of fine 
structures introduced into the quantum state. The scale at which these structures vary allows us to 
identify the minimal requirements for the precision of the complete interferometric sequence. Our 
results—especially the explanation of the underlying principle of the entangling method—may help 
to develop ultra-precise interferometers. 

PACS numbers: 37.25.-l-k, 03.75.-b,67.85.Hj 


Introduction —Generation of many-body entangled 
states for future applications in precise interferometry is 
a challenging task. With photons, it can be accomplished 
by creating a squeezed state of light which is fed into one 
of the interferometer’s arms [ij. Here squeezing refers to 
the reduction of the fluctuations of one of the quadratures 
of the electromagnetic field. The concept of squeezing 
can be—to some extent—transferred to the atomic do¬ 
main. However, since coherences between different atom 
number states are forbidden Q, a single-mode squeezed 
state of massive particles does not exist. Nevertheless, 
squeezing—in this context called spin-squeezing—can be 
generated by reducing the fluctuations of the atom num¬ 
ber difference between the two arms of the interferometer 
ii- This reduction is usually achieved using the repul¬ 
sive two-body interactions between the atoms [^. 

In the pioneering work Q , spin-squeezing was observed 
in a collection of ultra-cold rubidium atoms trapped in 
a double-well potential. The following experiments em¬ 
ployed a technique called one-axis twisting (OAT) to 
squeeze the fluctuations between two hyper-fine states of 
ultra-cold atoms [ffl- In this approach, first discussed 
in Q, particles in a coherent two-mode state 131 are 
dynamically entangled due to collisions. Spin-squeezing 
emerges immediately as the interactions start to operate. 
At some moment, however, spin-squeezing is lost, yet this 
does not mean the state is useless for precise interferom¬ 
etry. To contrary, the amount of useful entanglement 
continues to grow 1^, reaching the maximal attainable 
value [l^ . This effect—loss of squeezing accompanied by 
the increase of the entanglement useful for interferometry 
witnessed by the Fisher information—has been recently 
observed in the experiment [l^ . 


Interferometry with matter-waves seems particularly 
promising for future realistic metrological applications 
[l7| . This is because massive particles can strongly inter¬ 
act with the external field imprinting the phase, giving 
high interferometric signal. Moreover, two-body colli¬ 
sions of atoms allow for creation of many-body entangled 
states without the mediating element, such as a crystal 


used for the parametric down conversion of light 


In this work, we identify the underlying mechanism re¬ 
sponsible for high interferometric efficiency of the OAT 
states. We consider two emblematic transformations— 
the Mach-Zehnder interferometer (MZI) and a general¬ 
ized beam-splitter (BS)—both accessible with current ex¬ 
perimental techniques 0i, E 0, nil- We show that 
these interferometers, combined with the measurement of 
the number of particles in each output port, provide an 
ultra-high phase sensitivity. We also discuss the influence 
of the detection imperfections. Those come from the lim¬ 
ited control over the interaction strength or the duration 
of the OAT procedure as well as the imprecise measure¬ 
ment of the number of particles. Our results contribute 
to the understanding of the general features of the quan¬ 
tum nonlinear rotator (^ . They may also find applica¬ 
tion to the developing field of quantum-enhanced inter¬ 
ferometry operating on cold-atom ensembles E 231 and 
Bose-Einstein condensates. The experimental achieve¬ 
ment of [ifij l underlines the practical relevance of the 
entangled—but not spin-squeezed—states. The content 
of our work emphasizes the practical significance of even 
further-correlated OAT states for ultra-precise interfer¬ 
ometry. Other authors studied the role of decoherence in 
the OAT procedure (see for example [illii). We focus 
on the implementation of the OAT states and empha¬ 
size the role of the precise detection in maintaining high 
sensitivity. 


Twisted states and quantum interferometry —Usually, 
an interferometer consists of two arms through which ei¬ 
ther a pulse of light- or a matter-wave moves. During this 
propagation, a relative phase 6 is imprinted between the 
arms which later intersect to give an interference signal. 
From this output, the value of 9 is deduced. 


The OAT is a method of preparation of an entan¬ 
gled state, which is injected into an interferometer and it 
works as follows. Take a pure Bose-Einstein condensate 
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in an equal superposition of modes a and &, i.e., 


|V>) = 


1 

^/w. 


(a) + b^\ 

\ ^ ) 


N 

|o), 


( 1 ) 


where and w are the corresponding creation opera¬ 
tors. The goal is to transform this separable state of N 
particles into an entangled one. This is achieved in the 
absence of tunneling between the modes, when the evolu¬ 
tion is triggered by the two-body interactions. We denote 
by X~^ the associated timescale and have the two-mode 
Hamiltonian 


H = hxJl (2) 

Here Jz is one member of a triad of angular momentum 
operators Jx = \{a)b fetd), Jy = ^(d^d — b^a), and 
Jz = \{a)a — wb). The evolution of the initial coherent 
state m with the Hamiltonian ([5]) gives 

(3) 

where a = xt- The effect of interactions can be pictured 
using the mode occupation states |m), which denote ^ -I- 
m particles in mode a and ^ in ^ with m G [“ y; y]- 
It this basis, the state ([3|) reads 

JV 

IV'a)= E (4) 


with Cm = . The action of the Hamiltonian 

@ imprints the a-dependent phase which is quadratic 
in m. For most values of a, this phase jumps rapidly 
between the states \m) and |m-|-1). This is an important 
property of the OAT, and it is useful for the analytical 
calculations. 

The OAT procedure transforms the coherent state o 
into the entangled one ([3]) and is usually considered in 
the relation to the spin-squeezing parameter ii 


£2 _ Ar 


( 5 ) 



FIG. 1. (color online) The QFI from Eq. ([8| optimized over 
the direction of n (dashed black line) compared to the QFI for 
the BS (solid blue line) and for the MZI (solid red line). The 
inverse of calculated using Eq. ([5]) optimized over the di¬ 
rection of ni is shown with a green dotted line. All quantities 
are drawn for N = 100 particles. 


^± 2 ima^ Apart from a ^ kir with fc G Z, these phase 
factors quickly oscillate. As a result, the denominator of 
Eq. ([5]) is strongly suppressed. This means that within 
one period of evolution the state can be significantly spin- 
squeezed only for small a (i.e., short times), as observed 
in the experiment 1^ and in agreement with the numer¬ 
ically calculated which optimized over the direction 
of fii at each a is drawn in Fig. [T] with dotted green line. 
Analytical estimates show that indeed the OAT gener¬ 
ates spin-squeezed states only for a < TV “2 [1^, when 
the phase variations are smooth on the m —>■ to -I- 1 in¬ 
crement 27|. This means that for large N the presence 


of spin-squeezing does not last long. 

Although > 1 at later times, the state ([3]) is still 
useful for SSN interferometry IJ, l2l|. To demonstrate 
this, note that an interferometer imprints 9 through a 
unitary transformation 


H„(0) = e-”'^". (6) 


Each Jm = J ■ Hi is a scalar product of the vector 

of angular momentum operators J = [Jxt Jy, Jz] and a 
unit vector Hi {Hi, H 2 , and H 3 are orthogonal), while 
AJ„j = Jnj — (Jni)- When < 1, the state is called 
spin-squeezed, and it follows that it is particle-entangled 
and useful for sub shot-noise (SSN) metrology 14 l26l|. 
This is because the phase sensitivity can be linked to the 
spin-squeezing by the formula A9 oc ^/y/N. Note that 
< 1 can happen only when the denominator of Eq. dSD 
is large. Since the angular momentum operators couple 
the neighbouring states |to) and |to -I- 1), then, accord¬ 
ing to Eq. o, their average value contains phase factors 


Similarly to Eq. ([5]), = J ■ H while different direc¬ 

tions of H refer do distinct interferometric transforma¬ 
tions. When H is oriented along the x axis, thus = Jx, 
the mode mixing transformation represents a generalized 
BS, while for Jn = Jy, it stands for the MZL Finally, 
for Jn = Jz, the interferometer only imprints the phase 
9 between the modes without mixing them. Any other 
direction of H can be expressed as a sequence of these 
transformations. 

If 9 is estimated through a series oi 1 experiments, 
the precision of estimation, according to the Cramer-Rao 
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Lower Bound [1^, cannot be better than 

A9 = ^^=. (7) 

Here, the Fq is the quantum Fisher information (QFI) 
[i^ . which tells how much information about 9 can be 
encoded in a state g through a transformation ([6|). For 
a pure state, such as in Eq. 0, it is given by four times 
the variance of the generator J„ calculated in this state, 

FQ=4((AJ„f). (8) 


When Fq = N, the lower bound for the sensitivity is 
the shot-noise limit (SNL) A9 = According to 

Eq. 0, when the QFI is large, higher sensitivity can be 
reached, but Fq > A is only possible—allowing for SSN 
sensitivity—when the initial state is particle-entangled 
0. The QFI calculated for the state 0 is drawn with 
a dashed black line in Fig. [T] for each a maximized over 
the direction n. It has a characteristic shape—it grows 
from the value of Fq = N for the initially uncorrelated 
state o to reach a wide plateau approximately at a mo¬ 
ment when the spin-squeezing vanishes. Finally, through 
a narrow peak it strictly reaches the Heisenberg limit 
Fq = at a = 5 and then symmetrically drops to 
Fq = N at a = tt [iJ| . This plot shows that for the state 
0, Fq ^ N for all a. 

To better understand why the QFI stays at a constant 
and very high level for most of the time, we focus on two 
types of interferometric transformations such that can be 
implemented in the experiment, i.e., the BS and the MZI 
(rotations around the x and y axis, respectively). When 
a > As—so that the phase in Eq. O qmckly jumps 
between the neighboring states—we have |30| 


Fq=4((AJ,/,)2):^2 
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(9a) 

(9b) 


where /Sm = \J (A/2 + m + l){N/2 — m), which for large 
A changes slowly as compared to Cm’s and can be ap¬ 
proximated by Y' Moreover, when a is not too close 
to I-, the quickly oscillating phase “kills” the line (I9bl) . 
giving Fq = iA^ both for the BS and the MZI. This 
happens for most values of a, which explains the pres¬ 
ence of the wide plateau in Fig. [TJ Only when a ~ 
line (HEl) starts to contribute since the phase factor is 
then ^ g*27rm _ ppj. |.]^g contribution of the 

line (l9bl) is positive, so it will increase the value of the 
QFI up to the Heisenberg limit at a = However, the 
MZI due to the opposite sign in line (I9bll will experience 
a major drop of the QFI. These analytical estimates fully 
agree with the numerical results shown in Fig. [TJ 


Equipped with this knowledge, we show that the es¬ 
timation of 9 from the measurement of the population 
imbalance between the two output ports of either the 
MZI or the BS interferometer exploits the entanglement 
of the state 0. To this end, we calculate the related 
Fisher information (FI) given by the formula 



where the probability p{m\9) is the projection of the out¬ 
put state of the interferometer onto the state |m), i.e., 

p{m\9) = \{m\UniO)\ipa)\'^ ■ (H) 


The FI measures the amount of information about 9 con¬ 
tained in the probability p{m\9), and the corresponding 
CRLB, in analogy to Eq. o, reads 


A9 = 


1 1 


( 12 ) 


Since the QFI sets the maximal attainable sensitivity for 
any measurement, F{9) ^ Fq holds [2^. For both BS 
and MZI transformations—chosen by aligning the direc¬ 
tion of the interferometric rotation n by either the x or 
y axis in Eq. m and plugging the resulting probability 
into the Eq. cni) — we calculate the FI as a function of 9 
for each a. For instance, for the BS, the FI is 


2 , 

m sin((pm-t-i ^^m) II T (13a) 


l\' 

/5m —1 Sin((^m ^m —l) | Cm — 1 \ J 


(13b) 


In this expression. Cm = {'m\Ux{9)\4>a) and (pm is its 
phase. On the plateau, when this phase quickly jumps 
between the neighbouring states, the contribution to the 
FI resulting from the multiplication of the line (I13ap by 
(I13b|) will be negligible. Moreover, the squared sine func¬ 
tions can be approximated by its average value 1 /2, which 
leads to 


F{9) 


N_ 
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/5m I 


m— — 
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(14) 


This expression closely resembles Eq. ((9aJl and differs 
only by a constant factor of two. When 0 ~ 0, so that 
Cm remains peaked around m = 0, we can approximate 
/5^ ~ A^/4 which gives 


N 



(15) 


Note that for the MZI, in Eq. (HI, the sine functions are 
replaced with the cosine and the minus sign between lines 
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(I13a|) and (I13b|) is replaced by a plus. This however has 
no impact on the above discussion giving the same value 
of the FI as in Eq. m- Therefore, for both interferom¬ 
eters, the estimation from the population imbalance is 
very efficient and gives the Heisenberg scaling for a wide 
range of a’s. The outcome of our analytical considera¬ 
tions fully agrees with the numerical calculations of the 
FI directly from Eq. (Unj. 

The result ([T5|l is obtained for the ideal case, in absence 
of detection imperfections. We now make a final step and 
analyze the influence of such limitations. 

Influence of detection imperfections —We consider two 
sources of imperfections which might have a dominating 
impact on the performance of the OAT-based interfer¬ 
ometer. The first one results from an inadequate control 
over the value of a. This might be a consequence of the 
limited knowledge about the strength of the two-body 
interactions (for instance resulting from the fluctuating 
magnetic field (3ll - [33 |l or the lack of precise control over 
the duration of the OAT procedure. We account for this 
effect by assuming that a, rather than fixed, fluctuates 
with a Gaussian probability with uncertainty 6 a. In con¬ 
sequence, we replace the pure state from Eq. m with a 
mixture 



FIG. 2. (color online) The FI with finite atom detection res¬ 
olution (7, normalized to the FI with perfect accuracy. Three 
different solid lines correspond to different values of 0 ~ 0 
(i.e., 6 = 0.027r, O.OTtt and O.ldyr). The dashed black line is 
the best fit giving a function 0.095a~^. The horizontal dashed 
grey lines indicate the SNL for N = 100, 300 and 1000. 

F{a) calculated at the plateau for different values of 9, 
see Fig. [2j We observe that the value of FI has a universal 
behavior 


Qa — 


J daP{a 


afifa) {'Ip&l , 


(16) 


1 

F(a)cx —X-. 
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(19) 


with P{a) oc exp [—d^/2(da)^]. The FI from Eq. (fTUl) 
is now based on the probability distribution p{m\9) = 


Tr 


|m) (m| Uni9)g^Uli0) 


It is reasonable to assume 


that da < 1, otherwise it would be not possible to resolve 
even the wide plateau of Fig. [T] In this limit, the result 
of Eq. (ITKI) is replaced by 


^ 1 

~ ^ ^/l + 2N{6a^' 


(17) 


If the control over a is high so that da < N~ 5, the scaling 
proportional to is retained. On the other hand, when 
da > A “2 (or even is of the order of unity), the scaling 
is A oc As, which still is still much better than the SNL. 

Another factor which might reduce the sensitivity is 
the finite resolution of the atom detection a . To 

account for this, we replace the probability (HU with 


p{m\0) 


jV 

2 

R{m — m!)p{m'\0), 



(18) 


where we take R{m) oc exp[—m^/2cr^] as the resolution 
function. The finite resolution may have a much more 
profound effect on the FI than the uncertainty in de¬ 
termination of a. This is because on the plateau, its 
characteristic rapid phase variations summed even over 
a small increment of m' oscillate down to zero. This in¬ 
tuitive picture is confirmed by the numerically obtained 


This cr-dependent scaling function does not depend on 
A because it is rather related to the fast variation of the 
phase on the m' ^ m' + 1 increment in Eq. (fTSl) . For 
given A, this expression allows for a simple estimate of a 
minimal required resolution to beat the SNL. In general, 
^ is necessary to retain the SSN sensitivity, see 

the horizontal dashed lines in Fig. [51 

The above result is an example of a more general fea¬ 
ture of entangled states. From the metrological point of 
view, such states are highly susceptible to the parameter¬ 
imprinting transformations [s^. This property is wit¬ 
nessed by the fidelity between two neighboring prob¬ 
abilities p(to|0) aiidp{m\9 + S0) [l^. Generally speaking, 
probabilities which have fine structures are more prone 
to changes when 9 is varied. In our case, those fine struc¬ 
tures can be found in the rapid variations of the phase, 
which are governed by the value of a. It is now clear that 
on the plateau, when a is high, so that the phase changes 
abruptly, a sophisticated detection technique is required 
to grasp the quantumness of the state IV'a)- 

Conclusions —We have shown that the rapid phase 
variations of the one-axis-twisted state is the underlying 
mechanism responsible for its high interferometric effi¬ 
ciency. We have focused on two experimentally acces¬ 
sible transformations. The measurement of the popula¬ 
tion imbalance between the output ports provides a high 
sensitivity, which might be altered by the experimental 
imperfections. The lack of knowledge about the precise 
value of a is not as severe, as compared to the inaccu- 














5 


rate atom number detection. The presented discussion, 
in particular the explanation of the essential principle of 
the one-axis-twisting method, might find application in 
other areas of ultra-precise metrology. 
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